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Abstract 

Wilson loop averages are evaluated for large contours and in the large- A'c limit by 
means of minimal surfaces. This allows the study of the quark-antiquark gauge invariant 
Green function through its dependence on Wilson loops. A covariant bound state equation 
is derived which in the center-of-mass frame and at equal-times takes the form of a Breit- 
Salpeter type equation. The interaction potentials reduce in the static case to a confining 
linear vector potential. For moving quarks, flux tube like contributions are present. The 
nonrelativistic limit is considered. 
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The Wilson loop approach is based on the use of the gluon field path-ordered phase 
factor along a line C joining a point x to a point y: 

ry 

-ig / dzf^AJz) 

U{Cy,,y,x)^U{y,x)=Pe A ^ (1) 

Equations satisfied by phase factors were obtained and analyzed by Mandelstam j2] and 
Nambu 

The Wilson loop is defined as the trace in color space of the phase factor on a closed 
contour C: 



*(C) = ^trPe ■'Jc (2) 



-ig <j> dx^A^{x 

Its vacuum expectation value is denoted W{C): 

W{C) = mC)). (3) 

Loop equations were obtained by Polyakov and Makeenko and Migdal IHl El El • The 
Wilson loop essentially satisfies two types of equation, which are equivalent to the QCD 
equations of motion: 

1) The Bianchi identity. 

2) The loop equations (or Makeenko-Migdal equations). Those actually represent an 
infinite chain of coupled equations. 

A third property, factorization, is obtained in the large- Ai'c limit for two disjoint con- 
tours: 

W{CuC2) = W{C,)W{C2). (4) 

Considerable simplification occurs in the large- A*";, limit of the theory 1^1. In that 
limit, for large contours, i.e., at large distances, nonperturbative asymptotic solutions to 
the Wilson loops are represented by the minimal surfaces having as supports the loop 
contours ^ Among various types of surfaces, minimal surfaces are the only ones 
that satisfy the Bianchi identity; actually the latter becomes identical to the defining 
equation of minimal surfaces Furthermore, if at large distances one ignores short- 
distance perturbative solutions, then the unrenormalized coupling constant which appears 
in the loop equations can be adjusted in relation with the string tension ^01 and therefore 
minimal surfaces define on their own independent solutions. That property fails, however, 
when short- distance perturbative solutions are taken into account, since in this case the 
unrenormalized couphng constant is adjusted with respect to the perturbative regime. 
Minimal surfaces then become only large-distance asymptotic solutions. Nevertheless, 
if one is interested only in the large-distance behavior of the theory, saturation of the 
Wilson loop averages by minimal surfaces provides a correct description of the theory in 
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this regime. In that case, the Wilson loop average can be represented by the following 
functional of the contour C: 

=e-'^^(^), (5) 

where a is the string tension and A{C) the minimal area with contour C . 

Minimal surfaces also appear as natural solutions to the Wilson loop averages in two- 
dimensional gauge theories [TT] . 

To deal with the quarkonium bound state problem, we start with the two-particle 
gauge invariant Green function for quarks qi and q2 with different flavors and with masses 
mi and m^. 

G{XI,X2]X[,X'2) = {ip2ix2)U{x2, Xi)ipi{Xi)ipi{x[)U{x[, X2)^2ix2)) A,gug2- (6) 

Here U{x2, Xi) is the phase factor taken along the straight-line X1X2 (and similarly for 
U{x[,x'2)). Integrating in the large- iVc limit with respect to the quark fields, one obtains: 

G{xi,X2;x[,x'2) = -{tlU{x2,Xi)Si{xi,x[)U{x[,X2)S2{x'2,X2))A, (7) 

where 5*1 and 5*2 are the quark and antiquark propagators in the presence of the external 
gluon field and tr designates the trace with respect to the color group. 

The Green function G satisfies the following equation with respect to the Dirac oper- 
ator of particle 1 acting on xi. 

ih-d^xi) - mi)G{xi,X2;x[,X2) = -z(tr f/(x2, - x[)U{x[,x'2)S2{x2,X2))a 

JO dx"{a) 

where the segment X1X2 has been parametrized with the parameter a as x{a) = (1 — 
+ ax2', furthermore, the operator d/6x°' does not act on the explicit boundary point 
Xi of the segment, this contribution having been cancelled by the contribution of the 
gluon field A coming from the quark propagator Si. A similar equation also holds with 
the Dirac operator of particle 2 acting on X2- The operation SU{x2,xi)/6x{a) introduces 
an insertion of the field strength F at the point x{a) of the straight-line X1X2 I21I11IH]- 

In order to make apparent the Wilson loop structure of the two-particle Green function, 
we adopt a representation for the quark propagator in the external gluon field based on an 
explicit use of the phase factor along straight lines fU]- Introducing the gauge covariant 
composite object S{x,x'), made of a free fermion propagator So{x — x') (without color 
group content) multiplied by the path-ordered phase factor U (x, x') [Eq. ([T])] taken along 
the straight segment x'x, 

S{x,x') = So{x — x')U{x,x'), (9) 
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one shows that the quark propagator S{x, x') in the external gluon field satisfies the 
following functional integral equation in terms of S: 

S{x, x') = S{x, x') - I d^x"S{x, x")7" dX X—^S{x'\ x'), (10) 



5x"(A) 

where the operator 6/6x"'{X) acts on the factor U of S, along the internal part of the 
segment x'x", with x' held fixed. A similar equation in which the roles of x and x' are 
interchanged also holds. Those equations lead to iteration series for 5* in which the gauge 
covariance property is maintained at each order of the iteration. 

Using the above representations for the quark propagators in Eq. ((2)) one obtains for 
the two-particle Green function a series expansion where each term contains a Wilson 
loop along a skew-polygon: 

oo 

G = E G^J^ (11) 

where Gij represents the contribution of the term of the series having [i — 1) points of 
integration between xi and x[ {i segments) and (j — 1) points of integration between X2 
and X2 (j segments). We designate by Cij the contour associated with the term Gij. A 
typical configuration for the contour of 6*4,3 is represented in Fig. [T] 




Each segment of the quark lines supports a free quark propagator and except for the 
first segments (or the last ones, depending on the representation that is used) the Wilson 
loop is submitted to one functional derivative on each such segment. One then uses for 
the averages of the Wilson loops appearing in the above series the representation with 
minimal surfaces [Eq. (0)]. 
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Representation (fTU]) for the quark propagator is also used in Eq. (jHl) and its partner 
satisfied by the two-particle Green function G. One obtains two compatible equations for 
G where the right-hand sides involve the series of the terms Gij of Eq. (fTTj) and their 
functional derivative along the segment xiX2- In order to obtain bound state equations, 
it is necessary to reconstruct in the right-hand sides the bound state poles contained in 
G In x-space, bound states are reached by taking the large separation time limit 

between the pair of points {xi,X2) and {x'^^.x^) [T^ . 

To produce a bound state pole, it is necessary that there be a coherent sum of con- 
tributions coming from each Gjj, since the latter, taken individually, do not have poles. 
It is evident from representation ^ that functional derivatives acting on Wilson loop 
averages give rise to functional derivatives of the corresponding minimal surfaces Ai^. 
One then separates the various contributions of the right-hand sides of Eq. (jSI) and of 
its partner into different categories having the property of irreducibility. The first cate- 
gory corresponds to terms in which the functional derivative along the segment X1X2 acts 
alone on a given A. In the second category, the derivative along X1X2 is accompanied by 
another derivative along one of the segments of the quark lines, the two acting on the 
same A. In the third category, the derivative along X1X2 and two other derivatives along 
segments of the quark lines act on the same A or on the product of two As, respecting 
the irreducibility property of kernels, in the sense that they are not parts of the series 
expansion of a factorized G, and so forth. 

Let us now consider the terms of the first category. One notices that the derivative 
along the segment X1X2 acts on areas Ai^j with contour Cjj- which are different from one 
term of the series to the other (the number of segments being different). To have a 
coherent sum of those contributions it is necessary to expand each such derivative around 
the derivative of the lowest-order contour Ci 1, represented in Fig. |21 




Figure 2: The lowest-order contour Ci^i and its minimal surface Ai^i. 

It is that term that can be factorized and can lead through the summation of the 
factored series to the reappearance of the Green function G and to its poles. The remaining 
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terms do not lead to pole terms. Similarly, in the second category of terms, the two 
derivative contributions containing the derivative along X1X2 should be expanded around 
the lowest-order contribution coming from the contours C2,i or Ci 2, and so forth. 

In general, the derivative of the areas along X1X2 depends among others on the slope 
of the areas in the orthogonal direction to X1X2, thus feeling the positions of the other 
points x'l and Xj. In order to obtain bound state equations depending solely on xi and 
X2 and not on the remote points x'l and X2, it is necessary to associate the slopes of the 
areas along the segment X1X2 with the quark momenta. Taking then the large separation 
time limit, one ends up with two covariant compatible bound state equations, in which 
the interaction kernels or potentials are given by various functional derivatives involving 
at least one derivative along the segment X1X2. 

The fact that one has two independent but compatible equations means that the 
relative time variable does not play here a dynamical role and could in principle be 
eliminated. That is done by taking the difference and sum of the two equations. One 
finally ends up with one single three-dimensional equation at equal-times in the center- 
of-mass frame, having the structure of a Breit-Salpeter type equation ^J^]. Keeping 
for the potentials the terms containing one functional derivative of the area Ai^i [Fig. |2], 
the equation takes the form ^2] 

Po - {hio + h2o) - 7io7i - 72072 ^2/.] ^(x) = 0, (12) 

where is a 4 x 4 matrix wave function of the relative coordinate x = X2 — X1 considered at 
equal times, Pq the center-of-mass total energy and hiQ and /i2o the quark and antiquark 
Dirac hamiltonians; the Dirac matrices of the quark (with index 1) act on ip from the left, 
while the Dirac matrices of the antiquark (with index 2) act on ip from the right. The 
potentials Ai and A2 are defined through the equations 

= " lo " 5^^' '^''^ = ^ r "5^' ^'^^ 

x{a') belonging to the segment xiX2- 

The time components of Ai and A2 add up in the wave equation. For their sum, one 
has the expression (in the cm. frame) 

E1E2 f / El 

Aio + A20 = err 



arcsm 




+ (e(Pio) - e(P2o))( 




Ei + E2 



(14) 
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Here, r = vx^, Ea = ym"^ + p^, a = 1, 2, with the quark masses, p the cm. momen- 
tum, p = (p2 — Pi)/2, L the cm. orbital angular momentum, and e(pio) and e(p2o) the 
energy sign operators of the free quark and the antiquark, respectively: 

e{pao) = ^, a =1,2. (15) 

The space components of Ai and A2 are orthogonal to x. The expression of Ai is (in 
the cm. frame): 

E1E2 ( r2 E1E2 t 



~ ''''e, + E2\2L^E, + E,^ 



-|^arcsm(-y- 




(16) 



Here, p* is the transverse part of p with respect to x: 



1 



P = p-x— x.p. (17) 
x"^ 

The expression of A2 is obtained from that of Ai by an interchange in the latter of the 
indices 1 and 2 and a change of sign of p*. 

For sectors of quantum numbers where = 0, the expressions of the potentials 
become: 

^10 + A20 = ^(e(pio) + <P2o)W, (18) 
A.^-^(i(i.. + i..)-ii..)p'-' 

= ^E^2^l^^' ^ ~ ^^^^ 

The potentials are generally momentum dependent operators and necessitate an ap- 
propriate ordering of terms. 

From the structure of the wave equation (fT^ and the expressions of the potentials, one 
deduces that the interaction is confining and of the vector type. However, compared to the 
conventional timelike vector potential, it has additional pieces of terms contributing to the 
orbital angular momentum dependent parts. A closer analysis of those terms shows that 
they can be interpreted as being originated from the moments of inertia of the segment 
X\X2 carrying a constant linear energy density equal to the string tension. The interaction 
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potentials are therefore provided by the energy-momentum vector of the segment joining 
the quark to the antiquark, in similarity with the color flux tube picture of confinement. 
An analogous equation had also been proposed by Olsson et al. on the basis of a model 
where the quarks are attached at the ends of a straight string or a color flux tube [TT)| ITT] . 
A similar conclusion had also been reached by Brambilla, Prosperi et al. on the basis of 
the analysis of the relativistic corrections to the nonrelativistic limit of the Wilson loop 

[niiniiini. 

For heavy quarks, one can expand equation ()12|) around the nonrelativistic limit. To 
order 1/c^, the hamiltonian becomes (in the cm. frame): 



2h(y f I In 1/1 1 \, 0,0 ^ / 1 1 \C7 

^ + crr — + — - o ^ + ^ (P ) + ^ ^ + ^ - 

2/i 7r ^mi 8 ^mf 1712^ 4 ^mf my r 

a / 1 1 1 \ 9 cr /L.Si L.SoN 

-7r(— + — )^ + 2^ + ^{-^ + 

or ^mf 7712 mim2^ Ir^rrii m 
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2a f 1 1 2(T/ 1 1 

[/i = 77111712/ {mi + 7712), Si and S2 are the spin operators of the quark and of the anti- 
quark.] Several remarks can be made at this stage. First, the hamiltonian is independent 
of spin-spin interactions. Second, purely orbital angular momentum dependent pieces 
(proportional to L^) are present, the origin of which is related to the contribution to the 
rotational motion of the system of the moments of inertia of the flux tube, represented by 
the straight segment joining the quark to the antiquark. Those terms were also obtained 
in Refs. ^Sj and fOl Ej- Third, two kinds of spin-orbit term are present. The first 
comes from the contribution of a conventional timelike vector interaction represented by 
the potential ar, which is the dominant part of the combination ^410 + ^20 [Eq. (fTH) ]. The 
second type comes from the contributions of the direct interactions of the momentum of 
the flux tube with the quarks, represented by the spacelike potentials Ai and A2 [Eqs. 
(fT^]. The latter terms induce negative signs to the spin-orbit couplings, in opposite di- 
rection to the former one, a feature which is also observed on phenomenological grounds 
for the large-distance effects in fine splitting. 

The relativistic corrections to the interquark potential arising from the Wilson loop 
were analyzed and evaluated in the literature by Eichten and Feinberg [23 ? Gromes [22]) 
Brambilla, Prosperi et al. 1191 1^ . Brambilla, Pineda, Soto and Vairo [22] • 

The Wilson loop approach was also used for the study of quarkonium systems by 
Dosch, Simonov et al. with the use of the stochastic vacuum model 

In conclusion, the saturation of the Wilson loop averages in the large-A'^c limit by 
minimal surfaces provides a systematic tool for investigating the large- distance dynamics 
of quarkonium systems. A complete study of those systems necessitates the incorporation 
of the short- distance contributions, taken into account by perturbation theory. 
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